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Abstract 

We analyze the impact of recent solar and atmospheric data in the determination of the neutrino 
oscillation parameters, taking into account that both the solar v e and the atmospheric may 
convert to a mixture of active and sterile neutrinos. We use the most recent global solar neutrino 
data, including the 1496-day Super-K neutrino data sample, and we investigate in detail the impact 
of the recent SNO neutral current, spectral and day/night data by performing also an analysis 
using only the charged current rate from SNO. We confirm the clear preference of the pure active 
LMA solution of the solar neutrino problem and obtain that the LOW, VAC, SMA and Just-So 2 
solutions are disfavored with a A% 2 = 9, 9, 23, 31, respectively. Furthermore, we find that the 
global solar data constrains the admixture of a sterile neutrino to be less than 44% at 99% C.L.. 
A pure sterile solution is ruled out with respect to the active one at 99.997% C.L.. By performing 
an improved fit of the atmospheric data, we also update the corresponding regions of oscillation 
parameters. We find that the recent atmospheric Super-K (1489-day) and MACRO data have a 
strong impact on constraining a sterile component in atmospheric oscillations: if the Uu is restricted 
to the atmospheric mass states only a sterile admixture of 16% is allowed at 99% C.L., while a 
bound of 35% is obtained in the unconstrained case. Pure sterile oscillations are disfavored with a 
Ax 2 = 34.6 compared to the pure active case. 

In the appendix we discuss the implications of the first 145.1 days of KamLAND data on the 
determination of the solar neutrino parameters. The inclusion of KamLAND enhances the rejection 
of non-LMA-MSW solutions by 13 units in A% 2 . The bound on the sterile neutrino fraction is 
practically unaffected in the boron-fixed case, while it improves from 61% to 51% in the boron-free 
case. 
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I. INTRODUCTION 



Apart from confirming, yet again, the long-standing solar neutrino problem P, Q, El 
the recent results from the Sudbury Neutrino Observatory (SNO) on neutral current 
(NC) events 0, 0| have given strong evidence that solar neutrinos convert mainly to an 
active neutrino flavor. In addition, valuable spectral and day/night information has been 
provided 0,0]. This adds to the already robust evidence that an extension of the Standard 
Model of particle physics is necessary in the lepton sector. Although certainly not yet 
unique, at least for the case of solar neutrinos, which can be accounted well by spin-flavor 
precession B i or „on -standard neutrino matter interactions [10J, the most popular joint 
explanation of solar and atmospheric experiments is provided by the neutrino oscillations 
hypothesis, with neutrino mass-squared differences of the order of Am^ < 10~ 4 eV 2 and 
AmJ TM ~ 3 x 10" 3 eV 2 , respectively. 

In the wake of the recent SNO NC results we have re-analyzed the global status of 
current neutrino oscillation data including these and the remaining solar data Q 3 E fl 
as well as the current atmospheric [JJJ, LL2J samples, including the 1489 days Super- 
Kamiokande data ^| and the most recent MACRO data Motivated by the stringent 
limits from reactor experiments we adopt an effective two-neutrino approach in which 
solar and atmospheric analyses decouple. However our effective two-neutrino approach is 
generalized in the sense that it takes into account that a li ght sterile neutrino 0, 0, 
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advocated to account for the LSND anomaly |22|, may take part in both 
solar and atmospheric conversions. The natural setting for such a light sterile neutrino 
is provided by four-neutrino models. In this paper we will determine the constraints on 
oscillation parameters in this generalized scenario following from solar and atmospheric data 
separately. Such separate analyses are necessary ingredients towards a combined analysis 
of all current oscillation data, including solar, atmospheric, negative short-baseline data 



and the LSND experiment [23, |24j. As shown in Ref. 23] such separate analyses can be 
performed independently of the details of the four-neutrino mass scheme. 

Since the release of the latest SNO data in April 2002 a number of g lobal solar neutrino 



analyses in terms of active oscillations appeared 3 1 . Moreover, 

it has been shown by model-independent comparisons of the SNO CC rate with the SNO 
NC and Super-K rates that transitions of solar neutrinos into sterile neutrinos are strongly 
constrained by the recent data (see, e.g., Refs. 0,0 0, 0])- However, so-far no dedicated 
global analyses exist, where a participation of a sterile neutrino in the oscillations is fully 
taken into account 1 . Here we present a complete solar neutrino analysis including sterile 



1 In R ef. l32l| admixtures of a sterile neutrino to solar oscillations are considered. However, the authors of 
Ref. [22| are mainly interested in the determination of the solar neutrino fluxes and hence, their results 
are complementary to those obtained here. Some considerations of sterile solar neutrino oscillations can 
also be found in Ref. 1301. 
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neutrinos, determining the allowed ranges for the oscillation parameters 9 S0L and Amj OL , as 
well as for the parameter < rj s < 1 describing the active-sterile admixture. Furthermore, 
we investigate in detail the impact of the SNO neutral current, spectral and day/night data 
and compare with an analysis where we use only the charged current rate from SNO. 
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Concerning the atmospheric data, we perform an update of previous analyses 
adopting again the most general parameterization of atmospheric neutrino oscillations in 
the presence of sterile neutrino mixing, characterized by four parameters. We find that 
the recent 1489-day Super-Kamiokande data combined with the latest MACRO data lead 
to considerably stronger rejection against a sterile neutrino contribution to the oscillations 
than the previous 1289-day data sample. 

The plan of the paper is as follows. In Sec. Ill Al we set the general parametrization for 
solar oscillations in the presence of active-sterile mixing. In Sec. IIIBl we briefly describe the 
solar neutrino data and their analysis. In Sec. Ill CI we present the results of our analysis, 
aimed at studying the impact of recent solar data in the determination of the solar neutrino 
oscillation parameters, assuming, as mentioned, that the v e may convert to a mixture of 
active and sterile neutrinos. We give the regions of oscillation parameters for different 
allowed r] s values, display the global behavior of Axg 0L (Am;? 0L ) and Axs OL (0 SOL ) , calculated 
with respect to the favored active LMA solution, and evaluate the impact of the SNO NC, 
spectral and day/night data. Present solar data exhibit a higher degree of rejection against 
non-LMA and/or non-active oscillation solutions, which we quantify, giving also the absolute 
goodness of fit (GOF) of various oscillation solutions. Our solar neutrino results are briefly 
compared with those obtained in other recent analyses in Sec. Ill Dl In Sec. IIII Al we set 
our notations for atmospheric oscillations in the presence of active-sterile admixture. In 
Sec. IIII Bl we briefly describe the atmospheric neutrino data and their analysis. In Sec. IIII CI 
we describe our results for atmospheric oscillation parameters in an improved global fit of 
current atmospheric neutrino data. We quantify the impact both of our improved analysis 
as well as that of the recent data in rejecting against the sterile oscillation hypothesis. We 
update the corresponding regions of oscillation parameters and display the global behavior 
of Ax^ TM (Am^ TM ) and Ax^ TM (# ATM ). We compare the situation before-and-after the recent 
1489-day atmospheric Super-K data samples and give the present GOF of the oscillation 
hypothesis. In Sec. IIIID| we briefly compare our atmospheric neutrino results with those of 
other analyses. Finally, in Sec. IIVI we present our conclusions. 



II. SOLAR NEUTRINOS 

A. Active-sterile solar neutrino oscillations 

In the following we will analyze solar neutrino data in the general framework of mixed 
active-sterile neutrino oscillations. In this case the electron neutrino produced in the sun 
converts into a combination of an active non-electron neutrino u x (which again is a combi- 
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nation of and v T ) and a sterile neutrino v s : 



v e -> a/ 1 - rj s v x + y/rfs v s . (1) 

The parameter 7] s with < r] s < 1 describes the fraction of the sterile neutrino partici- 
pating in the solar oscillations. Therefore, the oscillation probabilities depend on the three 
parameters Am 2 0L , 8 S0L and r] s . The natural framework of light sterile neutrinos partici- 
pating in oscillations are four-neutrino mass schemes, proposed to account for the LSND 
result 22| in addition to solar and atmospheric neutrino oscillations. F or p revious studies 
of solar neutrino oscillation in a four-neutrino framework see Refs. 0, El Eland for an 
exact definition of the solar parameters and adopted approximations see Ref. |23j ]. 



B. Data and analysis 

As experimental data, we use the solar neutrino rates of the chlorine experiment Home- 
stake (21 (2.56 ± 0.16 ± 0.16 SNU), the most recent result of the gallium experiments 
SAGEyJ (70.8 tU ±|;I SNU) and GALLEX/GNO fl (70.8 ± 4.5 ± 3.8 SNU), as well 
as the 1496-days Super-Kamiokande data sample ij in the form of 44 bins (8 energy bins, 
6 of which are further divided into 7 zenith angle bins). In addition to this, we include the 
latest results from SNO presented in Refs. 0,Q|, in the form of 34 data bins (17 energy bins 
for each day and night period). Therefore, in our statistical analysis we use 3 + 44 + 34 = 81 
observables, which we fit in terms of the three parameters Am^, and r] s , with a Xsol 
of the form 

81 

xL(AmL, SOL , Vs) = £ (rt - Rf) ■ (<4 + tikfc 1 ■ (W - R?) ■ ( 2 ) 

In order to fully isolate the impact of the recent neutral current, spectral and day /night 
information of the SNO result, we also present an analysis which does not include such 
information. To this aim we use only the SNO events with energy higher than 6.75 MeV, for 
which the NC component is negligible 0. We sum these events to a single rate, combining 
with CI, Ga rates and full Super-K data, as described above. This procedure is analogous 
to the pre-SNO-NC situation, except that we take advantage of the enhanced statistics on 
the CC rate provided by the new data. We will refer to this analysis as SNO^cf analysis 
and it contains 48 data points. The comparison with the analysis including the complete 
SNO data published this year (SNO^^c) allows us to highlight the impact of the SNO NC, 
spectral and day/night information. 

For the solar neutrino fluxes we use the Standard Solar Model (SSM) flux [3^, including 
its standard 8 B flux prediction 2 . Motivated by the excellent agreement of the recent SNO 



We choose not to include the flux indicated by the recent Sn measurement of Ref. |3 



4 



NC result with the predictions of the Standard Solar Model, we prefer to adopt a boron-fixed 
analysis. However, for case of the LMA solution we explicitly illustrate the effect of this 
assumption by performing also a boron-free analysis, where we treat the solar 8 B flux as 
free parameter in the fit. For simplicity we neglect the hep and F neutrino fluxes, whose 
contribution to the present solar neutrino experiments is marginal, while for the pp. Be, B, 
pep, N and O fluxes we u 8 e the SSM value giveu iu Ref. 0, taking properly into account 
their theoretical uncertainties and cross-correlations in the calculation of the x 2 function. 

For the neutrino cross sections of Chlorine, SAGE, GALLEX/GNO and Super-K we 
assume the same as used in previous papers 37, 38, 3^, while for the CC and NC neutrino 



deuteron differentia, cross sections relevant taSNO we use the tables given in Q. The 
contribution of the cross-section uncertainties to the covariance matrix for the Chlorine and 
Gallium experiments is calculated as suggested in Ref. j32 |. For a given experiment (Chlorine 
or Gallium) we use full correlation of the error on the cross section for low-energy neutrino 
fluxes (pp, pep, Be, N and 0), but no correlation of the cross section error between the 
low-energy fluxes and the higher-energy 8 B flux. 

The neutrino survival probability P ee is extracted from the neutrino evolution operator 
U, which we factorize as a product of three factors U sun , U vac and U eart h corresponding to 
propagation in the Sun, vacuum, and Earth, respectively. The first and last factors include 



matter effects with the corresponding density profiles given in Refs. |35( and |41|. As a 
simplifying approximation, we assume that U sun depends only on the neutrino production 
point x , U vac only on the Sun-Earth distance L and U eart h depends only on the zenith-angle 
£ of the incoming neutrinos. Therefore in our calculations we neglect the small correlation 



between seasonal effects and day-night asymmetry [42]. For each value of the neutrino 
oscillation parameters Am^/E, # SOL and r] s we calculate the neutrino survival probability 
P ee by averaging over xq, L and (, properly accounting for all the interference terms between 
U sun , U vac and U ear th- 

Special care is taken in including all the theoretical and experimental errors and their 
cross-correlations in the calculation of the covariance matrix, for which we follow the descrip- 
tion of Ref. 0] (covariance approach). In particular, the errors associated to the Boron-flux 
shape, the energy-scale and the energy-resolution uncertainties of the Super-Kamiokande and 
SNO experiments are recalculated for each point in parameter space. 



C. Results and discussion 

In order to determine the expected event numbers for the various solar neutrino ex- 
periments we calculate the u e survival probability for each point in parameter space of 
(tan 2 9 S0L , Amj 0L , i] s ) and convolute it with the Standard Solar Model neutrino fluxes ji^ 
and the relevant neutrino cross sections. We have compared such expected event numbers 
with the data described above, taking into account the detector characteristics and appro- 
priate response functions. Using the above-mentioned Xsol we have performed a global fit 
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Figure 1: Allowed regions of tan 2 9 SOL and Am^ for rj s = (active oscillations), r/ s = 0.2 and 
rj s = 0.5. The lines indicate the regions determined by the SNOqq analysis (see definition in text), 
the shaded regions correspond to SNO^^c (see text). The confidence levels are 90%, 95%, 99% 
and 3a for 3 d.o.f.. 

of solar neutrino data, whose results we now summarize. 
Our global best-fit point occurs for the values 

tan 2 6 SOL = 0.46, Am 2 OL = 6.6 x 10~ 5 eV 2 (3) 

and correspond to r] s = 0. We obtain a xLin = 65.8 for 81 — 3 d.o.f., leading to the 
excellent goodness of fit of 84%. In Fig. ^we display the regions of solar neutrino oscillation 
parameters for 3 d.o.f. with respect to this global minimum, for the standard case of active 
oscillations, rj s = 0, as well as for r) s = 0.2 and tj s = 0.5. The first thing to notice is the 
impact of the SNO NC, spectral, and day/night data in improving the determination of 
the oscillation parameters: the shaded regions after their inclusion are much smaller than 
the hollow regions delimited by the corresponding SNOqcT confidence contours. Especially 
important is the full SNO C q^ c information for excluding maximal solar mixing in the LMA 
region and in closing the LMA region from above in Am 2 OL . Values of Am 2 0L > 10 -3 eV 2 
appear only at 3a. Previously solar data on its own could not close the LMA region, only 
the inclusion of data from reactor experiments ruled out the upper part of the LMA 
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Figure 2: Ax 2 0L as a function of Airij 01 and tan 2 9 SO l, for pure active (rj s = 0), pure sterile (rj s = 1) 
and mixed neutrino oscillations (r] s = 0.5). Upper and lower panels correspond to the SNOqq 6 and 
SNO cc ' N q samples denned in text. 

region j3^|. We obtain the following 3a ranges (1 d.o.f.): 

LMA: 0.26 < tan 2 # SOL < 0.85, 2.6 x 10~ 5 eV 2 < Am 2 OL < 3.3 x 10~ 4 eV 2 . (4) 

It is interesting to note that these 3a intervals are essentially unchanged if we minimize 
with respect to r] s or if we apply the constraint r] s = (pure active oscillations). In order to 
compare our allowed regions given in Fig. ^ with those of other groups, one has to take into 
account that we calculate the C.L. regions for the 3 d.o.f. tan 2 ^sql, Am s 2 0L and r] s . Therefore 
at a given C.L. our regions are larger than the usual regions for 2 d.o.f., because we also 
constrain the parameter rj s . 

Next we notice the enhanced discrimination against non-LMA solutions implied by the 
new data, apparent in Figs. El and El This shows that the first hints in favor of a 

globally preferred LMA oscillation solution which followed mainly from the flatness of the 
Super-K spectra, have now become a robust result, thanks to the additional data, to which 
SNO has contributed significantly 3 . One sees that, in contrast with the SNO^f situation, 

3 See also Ref. Q- 
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SNOgg 5 


sno sp,dn 


Region 


tan 2 9 S0L Am 2 S0L X 2 SO l GOF 


tan 2 (9 SOL Am 2 0L Xsol GOF 


Pure active (ry s = 0) 


LMA 

LOW 

VAC 

SMA 

Just-So 2 


0.46 7.2 x 10~ 5 40.9 69% 
0.83 4.8 x 10" 8 46.2 46% 
1.7 6.6 x 10- 10 45.0 51% 
1.1 x 10~ 3 5.0 x 10~ 6 57.8 11% 
1.0 5.5 x 10~ 12 59.6 9% 


0.46 6.6 x 10~ 5 65.8 86% 
0.66 7.9 x 10" 8 74.4 62% 
1.7 6.3 x 10~ 10 74.4 63% 
1.4 x 10~ 3 5.0 x 10" 6 89.3 20% 
1.0 5.5 x 10~ 12 96.8 8% 


Mixed fag = 0.5) 


LMA 

LOW 

VAC 

SMA 

Just-So 2 


0.46 7.6 x 10~ 5 45.4 50% 
0.91 3.5 x 10" 8 51.1 28% 
1.6 6.9 x 10" 10 49.4 34% 
3.6 x 10" 4 4.0 x 10~ 6 59.7 8% 
1.0 5.5 x 10~ 12 59.8 8% 


0.42 6.6 x 10~ 5 74.4 62% 
0.83 4.8 x 10~ 8 86.3 27% 
0.35 4.6 x 10" 10 81.3 41% 
4.4 x 10" 4 4.0 x 10~ 6 96.3 9% 
1.0 5.5 x 10- 12 97.0 8% 


Pure sterile (rj s = 1) 


LMA 

LOW 

VAC 

SMA 

Just-So 2 


0.44 1.6 x 10" 4 56.0 15% 

1.6 1.4 x 10~ 9 58.5 10% 

1.7 6.9 x lO" 10 56.1 15% 
3.5 x 10" 4 3.5 x 10~ 6 61.2 7% 

1.1 5.5 x 10~ 12 59.9 8% 


0.38 1.6 x 10~ 4 99.0 6% 
1.6 1.1 x 10" 9 101.6 4% 
0.33 4.6 x 10~ 10 89.1 21% 
3.6 x 10~ 4 3.5 x 10~ 6 99.4 6% 
1.0 5.5 x 10~ 12 97.2 8% 



Table I: Best-fit values of Am 2 0L and 9 S0L with the corresponding Xs 0L and GOF for pure active, 
pure sterile, and mixed neutrino oscillations. Results are given for the SNOq^ 6 (left column) and 
for the full SNO^q'^, analysis (right column). The relevant number of d.o.f. is 48 — 2 (81 — 2) for 
the SNO^ e (SNOcc'Jc) analysis. 

non-LMA solutions do not appear at 95% C.L.. However, the LOW and VAC solutions still 
appear at 99% C.L. for 3 d.o.f.. 

In order to concisely illustrate the above results we display in Fig. [21 the profiles of Ax 2 0L 
as a function of Am 2 0L (left) as well as tan 2 (9 SOL (right), by minimizing with respect to the 
undisplayed oscillation parameters, for the fixed values of r] s = 0, 0.5, 1. By comparing top 
and bottom panels in Fig. [21 one can clearly see the impact of the full SNO^'™ sample 
in leading to the relative worsening of all non-LMA solutions with respect to the preferred 
active LMA solution. 

The corresponding best-fit values for the various solutions of Am 2 0L and (9 S0L and the 
values of Xsol evaluated at the best-fit points are compiled in Tab. [H This table gives 
results for the three cases considered above: pure active, pure sterile and mixed neutrino 
oscillations, both for the SNO^ 6 an d the full SNOj^j™ analysis. To calculate the goodness 
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Figure 3: Ax 2 OL displayed as a function of r] s with respect to favored active LMA solution, for the 
SNOq^ g (left panel) and the SNO^'™ (right panel) analysis, as defined in text. 



of fit of the various solutions we evaluate in this table the x 2 f° r 48 — 2 (81 — 2) d.o.f. for the 
SNO^cf (SNO cc ' NC ) analysis defined previously. Note that we fix i] s at the three values 0, 
0.5 and 1. In the pure active case we find for LOW, VAC, SMA and Just-So 2 the following 
differences in \ 2 relative to the global best-fit point in LMA 

AXlow = 8-7, A X 2 AC = 8.6, A X 2 MA = 23.5, A X 2 ust _ So2 = 31.0. (5) 

Note that especially SMA and Just-So 2 are highly disfavored with respect to LMA. 

In addition to the scrutiny of the different neutrino oscillation solutions in the solar 
neutrino oscillation parameters Am 2 0L and # SO l, the present solar data can test the sterile 
neutrino oscillation hypothesis, characterized by the parameter rj s introduced above. The 
results can be presented in several equivalent ways. For example, rejection of sterile solar 
neutrino oscillations is already hinted by comparing the middle and right panels of Fig. ^ 
with the left one, corresponding to the pure active oscillation case: clearly the solutions 
deteriorate as T] s increases. Furthermore, the lines for r] s = 0.5 and r] s = 1 shown in 
Fig. 121 clearly show that sterile solutions are strongly disfavored with respect to pure active 
solutions. 

In order to summarize the above results we display in Fig. 01 the profile of Ax 2 0L as a 
function of < r] s < 1, irrespective of the detailed values of the solar neutrino oscillation 
parameters Am 2 0L and 8 S0L . This figure clearly illustrates the degree with which the solar 
neutrino data sample rejects the presence of a sterile component for each one of the possible 
solar neutrino oscillation solutions. The figure shows how the preferred LMA status survives 
in the presence of a small sterile component characterized by r] s (also seen in Figs. ^ and 
I2J). Further, one sees that the value i] s = is always preferred, so that increasing rj s leads 
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Figure 4: Best-fit active solar neutrino survival probabilities. 



to a deterioration of all oscillation solutions. Notice that there is a crossing between the 
LMA and VAC solutions, as a result of which the best pure sterile description lies in the 
vacuum regime. However, in the global analysis pure sterile oscillations with 77^ = 1 are 
highly disfavored. We find a x 2 -difference between pure active and sterile of A% 2 _ a = 33.2 if 
we restrict to the LMA solution, or Ax 2 _ a = 23.3 if we allow also for VAC. For 3 d.o.f. the 
A%g_ a = 23.3 implies that pure sterile oscillations are ruled out at 99.997% C.L. compared 
to the active case. 

For the LMA solution we have also performed an analysis without fixing the boron flux 
to its SSM prediction. In this case we treat the 8 B flux as a free parameter in the fit, and 
remove the error on this flux from the covariance matrix. From Fig. El one can see that the 
constraint on rj s is weaker in the boron-free case than in the boron-fixed one, since a small 
sterile component can now be partially compensated by increasing the total 8 B flux coming 
from the Sun. From the figure we obtain the bounds 

solar data: r] s < 0.44 (boron- fixed) , r] s < 0.61 (boron-free) (6) 

at 99% C.L. for 1 d.o.f.. In summary, we have found that, as long as the admixture of 
sterile neutrinos is acceptably small, the LMA is always the best of the oscillation solutions, 
establishing its robustness also in our generalized oscillation scheme. 

To round off our discussion of the solar neutrino fit update we present in Fig. 0] the 
u e survival probability versus energy E for the various solutions LMA, LOW and VAC, 
calculated as described above at the local x 2 -minima given in Tab. |1J Similar plots can be 
made for the case of sterile oscillations. 



D. Comparison with other groups 

Before turning to the atmospheric neutrino fits let us compare our solar neutrino results 
with those of other groups. Since the release of the latest SNO data in April 2002 sev- 
eral analyses have appeared. Taking into account the large amount of experimental input 
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Fogli et al. [21], hep-ph/0206162 vl, Table 1 


Barranco et al. [9j, hep-ph/0207326 v2, Table 1 


present analysis 


d.o.f. 


75-3 


46 


75-3 


49-4 


80-3 


49-2 


41-4 


81-3 


81-3 


81-2 


81-2 


best fit 


LMA solution 


tan 2 6 S0L 


0.34 


0.38 


0.39 


0.41 


0.45 


0.45 


0.40 


0.41 


0.42 


0.47 


0.46 


Am 2 OL [1(T 5 eV 2 ] 


5.0 


6.9 


5.6 


6.1 


5.8 


7.9 


5.4 


6.1 


5.8 


5.6 


6.6 




57.0 


43.5 


50.7 


40.6 


75.4 


33.0 


30.8 


65.2 


73.4 


68.0 


65.8 


GOF 


90% 


58% 


97% 


66% 


53% 


94% 


80% 


85% 


63% 


81% 


86% 


^XlOW, active 


10.7 


9.0 


9.2 


10.0 


9.6 


8.1 




12.4 


10.0 




8.7 


^XvAC, active 




10.0 


25.6 


15.5 


10.1 


14. 




9.7 


7.8 




8.6 


^■XsMA, active 




15.4 


57.3 


30.4 


25.6 


23. 




34.5 


23.5 




23.5 


^^LMA, sterile 












29. 










33.2 


^XLOW, sterile 






















35.9 


^XvAC, sterile 










26.0 












23.3 


^XSMA, sterile 










39.7 












33.6 



Table II: Comparison of solar neutrino analyses among different groups. We show the number of 
analyzed data points minus the fitted parameters, the best-fit values of tan 2 9 S0L and Am 2 OL for 
active oscillations and the corresponding x 2 - m i n i m a and GOF. Further we show the A% 2 with 
respect to the best fit LMA active solution for various other solutions (active, as well as sterile). 



data, variations in the analysis (such as the construction of the y 2 function or the treat- 
ment of theoretical errors) and the complexity of the codes involved it seems interesting 
to compare quantitatively the outcomes of different analyses. In Tab. |H]we have compiled 
some illustrative results of the solar neutrino analyses performed by the SNO and Super-K 
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26 



27 



28, 



29 



30 



collaborations 0, H], as well as theoretical ones 0, EE, 

Generally speaking, on statistical grounds, one expects the differences in the statistical 
treatment of the data to have little impact on the global best-fit parameter values, which lie 
in the LMA region for all analyses and are in good agreement. These differences typically 
become more visible as one compares absolute values of the x 2 > and/or as one departs from 
the best-fit region towards more disfavored solutions. Aware of this, we took special care to 
details such as the dependence of the theoretical errors on the oscillation parameters, which 
enter in the covariance matrix characterizing the Super-K and SNO electron recoil spectra. 
This way we obtain results which we consider reliable in the full oscillation parameter space. 

In the row labeled "d.o.f." we show the number of analyzed data points minus the fitted 
parameters in each analysis 4 . One can see from these numbers that various groups use 
different experimental input data, in particular the spectral and zenith angle information of 
Super-K and/or SNO is treated in different ways. Despite obvious differences in the analyses 
there is relatively good agreement on the best-fit LMA active oscillation parameters: the 
obtained best-fit values for tan 2 # S0L are in the range 0.34 — 0.47 and for Am s 2 OL they lie in 
the interval (5.0 - 7.9) x 10" 5 eV 2 . There is also good agreement on the allowed ranges 
of the oscillation parameters (not shown in the table). For example, the 3cr intervals given 



in Ref. 
in Ref. 



25j (0.24 < tan 2 9 S0L < 0.89 and 2.3 x 10' 
3 (tan 2 9 SOL < 0.84 and 2.3 x 10~ 5 eV 2 



eV 2 < Am 2 0L < 3.7 x 10" 



< Am 2 OL < 3.6 x 10' 



eV 2 



4 eV 2 ) and 
agree very 



well with the ranges given in Eq. However, even for the favored LMA solution, there 



3 



to 



are some differences in the GOF of the best-fit LMA solution, ranging from 53% 
97% |27J], due to differences in the construction of the x 2 function by different groups. 

There is remarkable agreement on the rejection of the LOW solution with respect to LMA 
with a Ax 2 ow, active ~ 10- O ur result for the vacuum solution Axvac, acti 
agreement with the values obtained in Refs. [ll EH EH Ell > whereas Refs. (2 
higher values. Our result for the SMA solution of A^sma, a ctim = 23.5 is in good agreemen ; 
with the values obtained in Refs. 



8.6 is in good 
3, E3 obtain 



25l|30|,|3l(; while Refs. [26, 29] and especially Ref. |27| 



obtain higher values. On the other hand in Ref. |l| SMA is somewhat less disfavored 5 . 

There had been so-far no dedicated global analysis of solar neutrino oscillations including 
the most recent SNO data for the case where sterile neutrinos take part in solar oscillations 
(?7 S 7^ 0). Model-independent considerations of transitions into sterile neutrinos can be found 
in Refs. 



27 



_J. Solar neutrino oscillations in the presence of active-sterile admixtures 
are also studied in Ref. [^, although in a different context. In the lower part of Tab. [H] 
we compare the partial results given in Refs. |0] and 0| for the pure sterile case (r) s — 1) 
with the corresponding values found in the present analysis. Although there are noticeable 
differences of the shown Ax 2 -values, there is agreement on the qualitative behavior. We 



4 Here we do not treat rj s as a free fit parameter, since we consider only the limiting cases rj s = and 1; 
this is the reason for the number 81 — 2 in the present analysis. 

5 Tracing back the reason for these and other differences in Tab. ITT1 goes beyond the scope of this work. 
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have also studied intermediate levels of sterile neutrino admixture, giving the corresponding 
regions of oscillation parameters and the full x 2 profiles relative to the favored active LMA 
solution (not shown in Tab. HTl see Figs. [HE] and EJ). 

We now turn to the analysis of the latest atmospheric data. As already mentioned in 
the introduction, separate analyses of solar and atmospheric data samples constitute the 
necessary ingredients towards a full combined study of all current oscillation data, including 
also the short-baseline data, as shown in 23, 24 



III. ATMOSPHERIC NEUTRINOS 



A. Active-sterile atmospheric neutrino oscillations 



In our analysis of atmospheric data we will make use of the hierarchy Amg OL <C Am^ TM 
and neglect the solar mass splitting. Further, in order to comply with the strong con- 
straints from reactor experiments 15[ we completely decouple the electron neutrino from 
atmospheric oscillations 6 . In the following we will consider atmospheric neutrino data in a 
generalized oscillation scheme in which a light sterile neutrino takes part in the oscillations. 
The setting for such scenarios are four-neutrino mass schemes 16, 17, 3|. In such schemes, 
besides the solar and atmospheric mass-splittings, there is also a large Am 2 motivated by 
the LSND experiment In contrast with the case of solar v e oscillations, the constraints 
on the z/^-content in atmospheric oscillations are not so stringent: in fact such constrains are 



provided by atmospheric data themselves 
oscillations in this general framework [21 



4611^ As a result to describe atmospheric neutrino 
we need two more parameters besides the 



standard 2- neutrino oscillation parameters 9 ATM and AmJ TM . We will use the parameters d M 



and d s already introduced in Ref. 23], and defined in such a way that (1 — d^) and (1 — d s ) 



correspond to the fractions of and v s participating in oscillations with Am^ TM , respec- 
tively. Hence, pure active atmospheric oscillations with AmJ TM are recovered in the limit 
d^ = and d s = 1. In four- neutrino models there is a mass scheme-dependent relationship 
between d s and the solar parameter r/ s . For details see Ref. j^J. 

We will also perform an analysis by imposing the constraint d^ = 0. In such "restricted" 
analysis the is completely constrained to the atmospheric mass states. Only in this limit 
the parameter d s has a similar interpretation as r] s introduced in the solar case. For d^ = 
we obtain that z/„ oscillates into a linear combination of v T and u s with Am 2 4TM : 



d a = 



d s v T + v/ 1 - d s v s 



(7) 



For a dedicated study of these issues see Ref. 
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B. Data and analysis 



For the atmospheric data analysis we use all the charged- current data from the Super- 
Kamiokande 0] and MACRO 14 1 experiments. The Super-Kamiokande data include the 



e-like and /i-like data samples of sub- and multi-GeV contained events (10 bins in zenith 
angle), as well as the stopping (5 angular bins) and through-going (10 angular bins) up- 
going muon data events. We do not use the information on v T appearance, multi-ring [i and 
neutral-current events since an efficient Monte-Carlo simulation of these data sample would 
require a more detailed knowledge of the Super Kamiokande experiment, and in particular 
of the way the neutral-current signal is extracted from the data. Such an information is 
presently not available to us. From MACRO we use the through-going muon sample divided 
in 10 angular bins \ . We did not include in our fit the results of other atmospheric neutrino 
experiments, such as the recent 5.9 kton-yr data from Soudan- 2 | 47l| . since at the moment 
the statistics is completely dominated by Super-Kamiokande Furthermore, some of 

the older experiments have no angular sensitivity, and thus can not be used to discriminate 
between active and sterile neutrino conversion, our main goal. 

Our statistical analysis of the atmospheric data is similar to that used in Ref . |22( , except 
that we now take advantage of the new Super-Kamiokande data and of the full ten-bin zenith- 
angle distribution for the contained events, rather than the five-bin distribution employed 
previously. Therefore, we have now 65 observables, which we fit in terms of the four relevant 
parameters Am^ TM , ATM , and d s : 

65 

xi TM (Ami TM , £ ATM , d„, d s ) = £ W - N\ h ) ■ (a 2 ex + a 2 th )r/ . (iVf - Nf) . (8) 



Concerning the theoretical Monte-Carlo, we improve the method presented in Ref. |3£ 
by properly taking into account the scattering angle between the incoming neutrino and the 
scattered leptqn directions. This was already the case for Sub-GeV contained events, however 
previously |22| we made the simplifying assumption of full neutrino-lepton collinearity in 
the calculation of the expected event numbers for the Multi-GeV contained and up-going-yU 
data samples. While this approximation is still justified for the stopping and thru-going 
muon samples, in the Multi-GeV sample the theoretically predicted value for down-coming 

is systematically higher if full collinearity is assumed. The reason for this is that the 
strong suppression observed in these bins cannot be completely ascribed to the oscillation 
of the down-coming neutrinos (which is small due to small travel distance). Because of the 
non-negligible neutrino-lepton scattering angle at these Multi-GeV energies there is a sizable 
contribution from up-going neutrinos (with a higher conversion probability due to the longer 
travel distance) to the down-coming leptons. However, this problem is less visible when the 
angular information of Multi-GeV events is included in a five angular bins presentation of 
the data, as previously assumed 
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C. Results and Discussion 



Folding together the atmospheric neutrino fluxes |49j|, our calculated neutrino survival 
probabilities including Earth matter effects with the profile of Ref. j^], and the relevant 
neutrino cross sections, we determine the expected event numbers for the various atmospheric 
neutrino observables, taking into account the appropriate detector response characteristics. 
Comparing with the data described in Sec. 1111 Bl we have performed a global fit of the 
atmospheric neutrino data using the above-discussed Xatm> following the same method used 
in Ref. [3^. We now summarize the main features of this fit. 

Our global best-fit point occurs at the parameter values 

sin 2 ATM = 0.49, Am 2 TM = 2.1 x 10" 3 eV 2 (best) (9) 

and d s = 0.92, = 0.04. We see that atmospheric data prefers a small sterile neutrino 
admixture. However, this effect is not statistically significant, also the pure active case 
(d s — 1, dfj, — 0) gives an excellent fit: the difference in \ 2 with respect to the best-fit point 
is only Ax 2 ct _ best = 3.3. For the pure active best-fit point we obtain 

sin 2 ATM = 0.5, Am 2 TM = 2.5 x 10~ 3 eV 2 (active) (10) 

with the 3<T ranges (1 d.o.f.) 

0.3 < sin 2 ATM < 0.7, 1.2 x 10~ 3 eV 2 < Am 2 TM < 4.8 x 10~ 3 eV 2 (active). (11) 

The determination of the parameters # AT m and Am ATM is summarized in Figs. El and |3 At 
a given C.L. we cut the xItm a ^ a A% 2 determined by 4 d.o.f. to obtain 4-dimensional volumes 
in the parameter space of (9 ATU , Am 2 TM , d^, d s ). In the upper panels we show sections of these 
volumes at values of d s = 1 and d^ = corresponding to the pure active case (left) and 
at the best-fit point (right). Again we observe that moving from pure active to the best-fit 
does not change the fit significantly. In the lower right panel we project away both d^ and 
d s , whereas in the lower left panel we fix d s = 0.5 and project away only d^. Comparing the 
regions resulting from 1489 days Super-K data (shaded regions) with those from the 1289 
days Super-K sample (hollow regions) we note that the new data leads to a slightly better 
determination of 8 ATM and Am 2 TM . However, more importantly, from the lower left panel we 
see, that the new data shows a stronger rejection against a sterile admixture: for d s = 0.5 
no allowed region appears at 3a for 4 d.o.f.. 

In Fig. |U1 we display the A% 2 with respect to the global best-fit point as a function of 
sin 2 # ATM (for both signs of 9 ATM ) and Am 2 TM , minimizing with respect to the other parameter, 
for different assumptions on the parameters d s and d^. In contrast to the solar case shown 
in Fig. |21 the atmospheric x 2 exhibits a beautiful quadratic behavior, reflecting the fact that 
the oscillation solution to the atmospheric neutrino problem is robust and unique. Notice 
again the significant worsening of the fit for the case of a sizable sterile neutrino admixture 
(see, e.g., the line corresponding to d s = 0.5). 
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Figure 5: Allowed regions of the parameters Sill 17 ATM 

and Am 2 ATM at 90%, 95%, 99% and 3a for 4 
d.o.f. and different assumptions on the parameters d s and (see text). The lines (shaded regions) 
correspond to 1289 (1489) days of Super-K data. 

In Fig.[7|we summarize the behavior of atmospheric \ 2 with respect to the parameters d s 
and dp. Indeed, the most striking result of the present improved analysis is the stronger re- 
jection we now obtain on the fraction of the sterile neutrino 1 — d s in atmospheric oscillations. 
Fig. El (b) clearly illustrates the degree with which the atmospheric neutrino data sample 
rejects the presence of a sterile component. On this basis one can place a model-independent 
atmospheric limit on the parameter d s , 

atmospheric data: 1 — d s < 0.35 (12) 

at 99% C.L. (1 d.o.f.). For the case of the restricted analysis, in which d^ = 0, we obtain 7 

d^ = 0: l-d s <0.16. (13) 

By comparing Eqs. (fT2"j) and (JT3J) we note the importance of taking into account the finite 
dp value in the analysis. 

Although there is no substantial change in the 99% C.L. bounds on 1 — d s due to the 
new Super-K data there is a huge effect for the case of sizable sterile neutrino admixtures, 



7 Note that in this case the C.L. regions should be defined with respect to the "restricted" best-fit point, 
which occurs for d s = 0.99, and not with respect to the global one. 
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Figure 6: Ax 2 TM as a function of sin 2 9 ATU (left) and Am 2 TM (right), using 1289 (upper) and 1489 
(lower) days of Super-K data, for the case of neutrino oscillations with arbitrary d s and tL, best-fit 
d s and d^, pure active and mixed active-sterile neutrino oscillations. 

d s < 0.5. In Tab. II I II we have compiled the best-fit values of sin 2 6^™, Am 2 TM , the x 2 
values and the GOF for the various atmospheric data samples for pure active and pure 
sterile oscillations. In the last column we give the difference in x 2 between active and sterile 
oscillation cases. Comparing these numbers for the 1289 and 1489 days Super-K samples 
we observe that all the new data except the Sub-GeV sample lead to a significant higher 
rejection against sterile oscillations. In combination with MACRO data the 1289 days Super- 
K gave a difference between pure sterile and active oscillations of A% 2 _ a = 17.8, whereas 
with the recent data we obtain 

A X 2 _ a = 34.6, (14) 

showing that pure sterile oscillations are highly disfavored with respect to the active ones 8 . 
Let us note that MACRO data give an important contribution to this effect: MACRO 

8 Here we should remark that this big improvement in constraining the sterile component - which is clearly 
visible also in the analyses presented by the Super-K collaboration itself - cannot be explained only by 
the improved statistics provided by the new data sample. The leading contribution comes instead from a 
change in the data themselves, which may indicate that some modification in the experimental efficiencies 
has been introduced. However, we have verified that such changes do not affect the theoretical prediction, 
since no difference between 1289 and 1489 days is visible in the Monte-Carlo of the Super-K collaboration. 
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Figure 7: (a) Allowed regions of the parameters d s and at 90%, 95%, 99% and 3a for 2 d.o.f.. 
The lines (shaded regions) correspond to 1289 (1489) days Super-K data. Further we show Axi TM 
as a function of d s (b) and d^ (c), minimized with respect to all other parameters. The subscript 
"R" refers to the restricted analysis with du = 0. 



alone disfavors the sterile oscillations already with Axg_ a = 9.0. These limits on the sterile 
admixture are significantly stronger than obtained previously |23| and play an important role 
in ruling out four-neutrino oscillation solutions in a combined global analysis of the LSND 
anomaly 2^|. Note, however, that in contrast with the case of d s , there is no substantial 



improvement in constraining the parameter due to the new data, as seen in Fig. [7| (c). 

In order to better appreciate the excellent quality of the neutrino oscillation description 
of the present atmospheric neutrino data sample we display in Fig.lHlthe zenith angle distri- 
bution of atmospheric neutrino events. Clearly, active neutrino oscillations describe the data 
very well indeed. In contrast, no oscillations can be visually spotted as being inconsistent. 
On the other hand conversions to sterile neutrinos lead to an excess of events for neutrinos 
crossing the core of the Earth, in all the data samples except sub-GeV. 



D. Comparison with other groups 



Let us briefly compare our atmospheric neutrino oscillation results with those of other 
groups. Apart from the analyses presented in Refs. j3, 0] there had been no other complete 
atmospheric neutrino analysis taking into account the most general structure of neutrino 
mixing in the presence of sterile neutrinos, characterized by four mixing parameters. In 
the analyses of Refs. 



13, 



50] the is restricted to the atmospheric mass states, which 
corresponds to the constraint d^ = in our parameterization. However, at the corresponding 
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Active (d/j, = 0, d s = 1) 


Sterile (d^ = 0, d s = 0) 




Data sample 


d.o.f. 


sin 2 9 Am 2 [eV 2 ] xLt GOF 


sin 2 6 Am 2 [eV 2 ] xL GOF 




Super-K-1289 days, improved 


SK Sub-GeV 
SK Multi-GeV 
SK Stop-// 
SK Thru-/i 
MACRO 


20-2 
20-2 
5-2 
10-2 
10-2 


0.50 2.1 • 10~ 3 14.9 67% 
0.50 1.8 • 10~ 3 6.4 99% 
0.50 4.2 • 10-3 1.2 76% 
0.29 6.3 • 10~ 3 5.3 73% 
0.50 2.4- 10~ 3 11.0 20% 


0.50 2.2 • 10~3 15.0 66% 
0.57 3.5 • 10-3 11.3 88% 
0.61 4.0 • 10~3 3.1 38% 
0.84 1.0 • 10- 2 7.8 45% 
0.96 9.4 • 10~3 20.0 1% 


0.1 
4.8 
1.9 
2.5 
9.0 


SK Contained 

Upgoing-/i 

SK+MACRO 


40-2 
25-2 
65-2 


0.50 2.0 • 10~ 3 21.4 99% 
0.50 3.3 • 10"3 19.2 69% 
0.50 2.7 • 10-3 41.7 98% 


0.54 3.0 • 10"3 26.9 91% 
0.72 4.2 • 10-3 32.8 8% 
0.56 2.8 • 10-3 59.4 60% 


5.5 
13.6 
17.8 


Super-K-1489 days 


SK Sub-GeV 
SK Multi-GeV 
SK Stop-// 
SK Thru-yU 
MACRO 


20-2 
20-2 
5-2 
10-2 
10-2 


0.50 1.9 • 10-3 9.0 96% 
0.50 1.3 • 10-3 10.2 93% 
0.50 2.8 • 10-3 1.5 69% 
0.50 3.5 • 10-3 6.3 61% 
0.50 2.4-10-3 11.0 20% 


0.51 2.0 • 10-3 9.0 96% 
0.57 3.5 • 10"3 18.4 43% 
0.75 2.8 • 10~3 6.9 8% 
0.84 6.7 • 10-3 16.0 4% 
0.96 9.4 • 10~3 20.0 1% 


0.0 
8.2 
5.4 
9.7 
9.0 


SK Contained 

Upgoing-/i 

SK+MACRO 


40-2 
25-2 
65-2 


0.50 1.5 • 10"3 19.3 99% 
0.50 3.0 • 10"3 18.9 71% 
0.50 2.5 • 10-3 40.2 99% 


0.54 3.0 • 10"3 28.1 88% 
0.75 3.2 • 10"3 40.8 1% 
0.61 2.7 • 10-3 74.9 15% 


8.8 
22.0 
34.6 



Table III: Atmospheric neutrino best-fit oscillation parameters for pure active and pure sterile 
oscillations for the various data samples. 



limiting cases our generalized analysis can be compared with the results of other works. Let 



5ol | is based on the 1289-days SK data sample (79.5 



us further note that the analysis of Ref. 

kton-yr) and in contrast to Refs. 13j,|50| we use also data from the MACRO experiment. 

First, we find very good agreement in the case of pure active oscillations: the agreement 
of our best-fit values given in Eq. (Jl"Uj) with those obtained by the Super-K collaboration 
(sin 2 # ATM = 0.5, Am 2 , = 2.5 x lO -3 eV 2 [13]) is excellent, with good agreement also with 
the results of Ref. |50j (sin 2 6 ATM = 0.41, Am 2 TM = 3 x lO -3 eV 2 ). Similarly, also the allowed 
ranges shown in the upper left panel of Fig. [5] compare very well with the ranges obtained 
in Refs. 0E3- 

This shows that the determination of the active atmospheric oscillation 
parameters is already rather stable with respect to variations in the analysis and inclusion of 
additional of data. Concerning admixtures of sterile neutrinos, we note that it is presently 
not possible to use information on u T appearance, multi-ring \i and neutral-current events 
outside the Super-K collaboration, because to simulate these data a detailed knowledge of the 
detector and the applied cuts is necessary. These classes of events should provide additional 
sensitivity towards rejecting a possible contribution of sterile neutrinos. Therefore, the fact 
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Figure 8: Zenith angle dependence of the /j-like data used in our fit. Further we show the predicted 
number of atmospheric neutrino events for best-fit, pure- active and pure-sterile oscillations and no 
oscillations. 



that the value of Ax;?_ a = 49.8 [13( between pure active and sterile oscillations obtained by 
the Super-K collaboration is higher than our value 34.6 given in Eq. ()14j) is understandable, 
since with the Super-K data accessible to us we have a reduced discrimination between active 
and sterile oscillations, based solely on the matter effects. 



IV. CONCLUSIONS 

Prompted by the recent data on solar and atmospheric neutrinos we have reanalyzed the 
global status of oscillation solutions, taking into account the that both the solar u e and the 
atmospheric may convert to a mixture of active and sterile neutrinos. In addition to the 
SNO neutral current, spectral and day/night (SNO C q°q) results we add the latest 1496-day 
solar and 1489-day atmospheric Super-K neutrino data samples. 

We have studied the impact of the recent solar data in the determination of the regions 
of oscillation parameters for different allowed r\ s values, displaying the global behavior of 
Axhol^^Iol) an d Axg 0L (t9 S0L ), calculated with respect to the favored active LMA solution. 
We have investigated in detail the impact of the full CI + Ga rates + Super-K spectra + 
the complete SNO cc ' n q data set, comparing with the situation when the this year's SNO 
data is left out. We confirm the clear preference for the LMA solution of the solar neutrino 
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problem and obtain that the LOW, VAC, SMA and Just-So 2 solutions are disfavored with 
a Ax 2 = 9, 9, 23 and 31, respectively, for the pure active case. In addition, we find that 
the global solar data sample constrains admixtures of a sterile neutrino to be smaller than 
44% at 99% C.L.. This bound is relaxed to 61% when the solar 8 B flux is treated as a free 
parameter. A pure sterile solution is ruled out with respect to the active one at 99.997% 
C.L.. For allowed sterile neutrino admixtures LMA is always the best of all the oscillation 
solutions. We remark, however, the existence of non-oscillation solutions [8l, 9j, 101. These 
will be crucially tested 

B0 at the up-coming KamLAND r eact or experiment & 
By performing an improved fit of the atmospheric data, we have also updated the cor- 
responding regions of oscillation parameters for the case where the atmospheric convert 
to a mixture of active and sterile neutrinos. We have displayed the global behavior of 
^xi T M(^ m ATM) an d AXatm(^atm) for different allowed values of the sterile neutrino admix- 
ture in the atmospheric channel. We have compared the situation before-and- after the recent 
1489-day atmospheric Super-K data samples and shown that the GOF of the oscillation hy- 
pothesis is excellent. We have found that the recent 1489-day atmospheric Super-K data 
strongly constrain a sterile component in atmospheric oscillations: if the is restricted to 
the atmospheric mass states only a sterile admixture of 16% is allowed at 99% C.L., while a 
bound of 35% is obtained in the unconstrained case. Pure sterile oscillations are disfavored 
with a Ax 2 = 34.6 compared to the active case. 
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Appendix: IMPACT OF THE KAMLAND RESULT 

In a recent paper the first results of the KamLAND collaboration became public j^. 
The KamLAND experiment is a reactor neutrino experiment whose detector is located at 
the Kamiokande site. Most of the v e flux incident at KamLAND comes from nuclear plants 
at distances of 80 — 350 km from the detector, making the average baseline of about 180 
kilometers, long enough to provide a sensitive probe of the LMA-MSW region. The target 
for the v e flux consists of a spherical transparent balloon filled with 1000 tons of non- 
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Figure 9: A% 2 as a function of Am 2 0L and tan 2 # SOL , for solar+Chooz, KamLAND alone, and 
solar+reactor (both Chooz and KamLAND) data. 



doped liquid scintillator, and the anti-neutrinos are detected via the inverse neutron /?-decay 
process u e + p — > e + + n. The KamLAND collaboration has for the first time measured the 
disappearance of neutrinos produced in a power reactor. They observe a strong evidence 
for the disappearance of neutrinos during their flight over such distances, giving the first 
terrestrial confirmation of the solar neutrino anomaly and also establishing the oscillation 
hypothesis with man-produced neutrinos. 

In this appendix (which does not appear in the published version of this paper) we analyze 
the implications of the first 145.1 days of KamLAND data on the determination of the solar 
neutrino parameters. The details of our theoretical Monte-Carlo and statistical analysis are 
given in Ref . ; in particular, the KamLAND x 2 -function is calculated assuming a Poisson 
distribution for the experimental data, as described in Sec. IV of that paper. 

The impact of the KamLAND result on 6 S0L and Am 2 for the case of pure active oscil- 



lations has already been discussed in detail in Ref. |54|, and the results are summarized here 
in Fig. 121 First of all, we note that non-LMA solutions, characterized by a very small value 
of either Am 2 0L (LOW, VAC, Just-So 2 ) or # S0L (SMA) are in disagreement with the evidence 
for v e disappearance observed in KamLAND. As a consequence, the relative quality of these 
solutions with respect to LMA is worsened by an extra A% 2 « 13 when the KamLAND data 
are also included in the analysis, so that Eq. (JSJ) is now replaced by: 

A x 2 ow = 21.7, A X 2 , AC = 21.6, A X 2 MA = 36.5, A x 2 ust . So2 = 44.0. (A.l) 

This led to the conclusion that LMA is presently the only allowed solution to the solar neu- 
trino problem. The global best-fit point occurs for pure-active oscillations, and is practically 
unaffected by the inclusion of KamLAND (cfr. Eq. (j3J)): 

tan 2 6 S0L = 0.46, Am 2 0L = 7.2 x 10~ 5 eV 2 . (A.2) 
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Figure 10: Allowed regions of tan 2 8 SOL and Am 2 0L for tj s = (active oscillations), rj s = 0.2 and 
rj s = 0.5. The lines correspond to the analysis of solar+Chooz data, while the shaded regions 
correspond to the combination of solar+reactor (both Chooz and KamLAND) data. Both lines 
and shaded regions refer to 90%, 95%, 99% and 99.73% confidence intervals for 3 d.o.f.. 



However, as can be seen from Fig. the LMA region is now split into two sub- regions, and 
from the left panel of Fig. El we see that a secondary minimum characterized by tan 2 $sol — 
0.42 and Am 2 0L = 1.5 x 10~ 4 eV 2 appears. The relative quality of this point with respect to 
the global best-fit point given in Eq. (fA~2jl is A% 2 = 5.7. The allowed 3a ranges for Am 2 0L 
and 9 S0L are (cfr. Eq. (JH)): 

f 5.4 x 10~ 5 eV 2 < Am 2 < 9.8 x 10" 5 eV 2 , 
0.29 < tan 2 6 S0L < 0.85, <^ A ~ ° t ' A.3 

\ 1.3 x 10- 4 eV 2 < Am 2 0L < 2.0 x 10~ 4 eV 2 . 1 ' 

Let us now consider the impact of KamLAND on the determination of rj s . From Fig. El 
we have a first indication that the bound on the fraction of sterile neutrino participating in 
solar neutrino oscillations is essentially unaffected by this experiment. This can be easily 
understood since KamLAND is only sensitive to the anti-neutrino survival probability P ee , 
and is therefore unable to discriminate between different oscillation channels. Taking into 
account that matter effects induced by the Earth mantle are practically negligible given the 
short distance traveled by the neutrinos in their flight between the source and the detector, 
it is straightforward to conclude that KamLAND is completely insensitive to r] s . 

In the left panel of Fig. [TT] we display the profile of Ax 2 0L and Ax 2 0L+KAM as functions of rj s , 
irrespective of the detailed values of the solar neutrino oscillation parameters Am 2 0L and 6* SO l- 
The dashed lines correspond to free S B flux, while for the solid lines the boron flux is fixed 
to its SSM prediction. As expected, for the boron-fixed case the inclusion of KamLAND 
is almost completely irrelevant for the determination of the sterile neutrino fraction, and 
for i] s < 0.8 no visible difference arises between the pre-KamLAND and post-KamLAND 
cases. However, looking at Fig. ^2 we see that for the boron-free case the inclusion of the 



23 




Figure 11: Left panel: Axg 0L and ^Xsol+kam displayed as functions of rj s , for both fixed and free 
S B flux. Right panel: la allowed regions for Amg 0L from the analysis of solar+Chooz data, for 
both fixed (green) and free (red) S B flux. The dark (light) gray horizontal belts correspond to the 
la (2a) allowed regions for Am| OL from the analysis of KamLAND data. 



KamLAND data is relevant. At 99% CL. we obtain the bounds (cfr. Eq. ©): 

solar+reactor data: r] s < 0.43 (boron-fixed), r] s < 0.51 (boron-free). (A. 4) 

To understand why fixing or not-fixing the 8 B flux leads to such a different behavior, we 
illustrate in the right panel of Fig. ^2 the dependence on 7] s of the la allowed range for 
Airig 0L from the analysis of solar data alone, both for the boron-fixed (green) and the boron- 
free (red) cases. These two bands should be confronted with the horizontal gray belts 9 , 
which correspond to the KamLAND la (dark gray) and 2a (light gray) allowed intervals for 

Concerning the boron- fixed case, we see that for rj s < 0.6 the Amg OL value preferred 
by solar data is almost insensitive to r] s , and in very good agreement with the KamLAND 
prediction. In this regime a non-zero value of r] s leads to a mild deficit in the expected 
number of NC and ES events, thus reducing the quality of the fit. When r] s exceeds ~ 0.6, 
this deficit become relevant, and in order to compensate it the best-fit point moves towards 
regions of the neutrino parameter space where the electron neutrino survival probability P ee 
is larger. From Fig.[T2]it is easy to understand that this correspond to larger values of Am;? 0L . 
In any case, as long as rj s < 0.9 the Am| 0L regions allowed at la by solar and KamLAND 
data still overlap, and this explains why in this regime there is no visible difference between 

A XS0L aIld A XsOL+KAM- 

The situation is different for the boron-free case. As for the previous case, an increase 
of T) s lead to a deficit in the expected number of NC and ES events, which can now be 
compensated by assuming a larger value of the 8 B flux. However, a larger boron flux results 



9 Note that the gray bands are perfectly horizontal since KamLAND alone is insensitive to rj s . 
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Figure 12: Solar neutrino survival probability P ee as a function of E u /Am^ 0L , calculated at the 
solar+reactor best-fit point given in Eq. (|A.2|) . 

in an excess of v e arriving at the detectors, so that now the experimental data favor regions 
of the parameter space where the electron neutrino survival probability P ee is smaller. This 
explains why the preferred value for Am 2 0L decreases as rj s increases. When rj s > 0.8 the 
favored value for the 8 B flux rapidly decreases, and Am 2 0L increases again. A consequence 
of this is that already for rj s > 0.4 the la regions for solar and for KamLAND data no 
longer overlap. This leads to a tension between the two data sets, which results in an excess 
of Axsol+kam over ^Xsol- ^ * s om y f° r veT Y large values of i] s that solar data return in 
agreement with KamLAND, and the two lines Ax 2 0L and Ax 2 0L+KAM merges again. 

In summary, in this appendix we have investigated the impact of the first 145.1 days of 
KamLAND data on the determination of the solar neutrino parameters. We have found 
that all non-LMA solution are now ruled out, and that the original LMA region is split into 
two relatively narrow islands around the values of Am 2 0L = 7.2 x 10~ 5 eV 2 (best fit point) 
and Am 2 0L = 1.5 x 10~ 4 eV 2 (local minimum). The bound on the sterile neutrino fraction 
rjs is practically unaffected in the boron-fixed case, but improves from 0.61 to 0.51 (at 99% 
CL.) in the boron- free case. 
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